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Abstract. The nature of fast magnetoacoustic and Alfven waves is investigated in a zero fj plasma. This gives an indication 
of wave propagation in the low p solar corona. It is found that for a two-dimensional null point, the fast wave is attracted to 
that point and the front of the wave slows down as it approaches the null point, causing the current density to accumulate there 
and rise rapidly. Ohmic dissipation will extract the energy in the wave at this point. This illustrates that null points play an 
important role in the rapid dissipation of fast magnetoacoustic waves and suggests the location where wave heating will occur 
in the corona. The Alfven wave behaves in a different manner in that the wave energy is dissipated along the separatrices. 
For Alfven waves that are decoupled from fast waves, the value of the plasma fj is unimportant. However, the phenomenon of 
dissipating the majority of the wave energy at a specific place is a feature of both wave types. 

Key words. MHD - Waves - Sun: corona 

1. Introduction 

The coronal heating problem remains a key unsolved problem 
in solar physics. While the coronal magnetic field is ultimately 
involved, there are many rival theories ranging from reconnec- 
tion models involving nanoflares and wave heating models in- 
volving phase mixing and resonant absorption. The reconnec- 

\ tion models either require the formation of many current sheets, 
due to random photospheric boundary motions that braid the 

^ magnetic field, or the collapse of null points. The wave heat- 

■ ing models rely on the generation of small length scale wave 

" motions in the corona. There is clear evidence from SOHO 
and TRACE observations of slow MHD waves (Berghmans 
& Clette, 1999; De Moortel et ai, 2000), fast MHD waves 
(Nakariakov et ai, 1999) and non-thermal line broadening due 
to Alfven waves (Harrison et ai, 2002). While there may be in- 
sufficient energy in these waves to heat the whole corona, their 
dissipation will contribute to the overall energy budget. This 
paper is concerned with the propagation of MHD waves in the 
neighbourhood of null points in a zero /3 plasma, giving an in- 
dication of how MHD waves behave in the low /3 plasma of the 
solar corona. 

The existence of null points is predicted on theoretical grounds 
and their importance lies in the fact that the Alfven speed is 
actually zero at that point. This important consequence will 
be utilised later Potential field extrapolations, using photo- 
spheric magnetograms to provide the field distribution on the 
lower boundary, suggest that there are always likely to be 
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null points in the corona. The number of such points does 
depend on the magnetic complexity of the photospheric flux 
distribution. Detailed investigations of the coronal magnetic 
field, using such potential field calculations, can be found in 
IBeveridge et al. (2002)| and |Brown and Priest (2001)1 

The fact that null points are a weakness in the magnetic field 
has been used to investigate how they collapse in response to 
boundary motions. This has been investigated by Craig and co- 
workers and other authors using analytical and numerical ap- 
proaches (Craig & McClymont, 1991; Hassam, 1992). The ba- 
sic aim is to move the field lines passing through the boundary 
in a particular manner in order to perturb the field. The result- 
ing field perturbations cause the null point to collapse to form 
a current sheet in which reconnection can release magnetic en- 
ergy. In these models the boundary motions move the field lines 
but do not return them to their original positions. Thus, the 
Poynting flux induced by the imposed motion (and then fix- 
ing the field after the motion is complete) accumulates at the 
resulting current sheet and provides the energy released in the 
reconnection. However, if the boundary motions are simply due 
to the passing of incoming waves through the boundary, then it 
is not clear that the null point need collapse and form a current 
sheet. If this is the case, then it is not clear if the energy in the 
wave, again due to the Poynting flux through the boundary, will 
dissipate or simply pass through one of the other boundaries. 

Waves in the neighbourhood of 2D null points have been in- 
vestigated by various authors. Bulanov and Syrovatskii (1980) 
provided a detailed discussion of the propagation of fast and 
Alfven waves using cylindrical symmetry. In their paper, har- 
monic fast waves are generated and these propagate towards the 
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null point. However, the assumed cylindrical symmetry means 
that the disturbances can only propagate either towards or away 
from the null point. Craig and Watson (1992) mainly consider 
the radial propagation of the m = mode (where m is the az- 
imuthal wavenumber) using a mixture of analytical and numer- 
ical solutions. In their investigation, the outer radial boundary is 
held fixed so that any outgoing waves will be reflected back to- 
wards the null point. This means that all the energy in the wave 
motions is contained within a fixed region. They show that the 
propagation of the m - wave towards the null point generates 
an exponentially large increase in the current density and that 
magnetic resistivity dissipates this current in a time related to 
log J]. Their initial disturbance is given as a function of radius. 
In this paper, we are interested in generating the disturbances 
at the boundary rather than internally. Craig and McClymont 
(1991, 1993) investigate the normal mode solutions for both 
m - and m Q modes with resistivity included. Again they 
emphasise that the current builds up as the inverse square of 
the radial distance from the null point. However, attention has 
been restricted to a circular reflecting boundary. This paper will 
investigate wave propagation in a zero /3 plasma in the neigh- 
bourhood of a simple 2D null point but for more general dis- 
turbances, more general boundary conditions and single wave 
pulses. This allows us to concentrate on the transient features 
that are not always apparent on using normal mode analysis. 
Galsgaard et al. (2003) looked at a particular type of wave dis- 
turbance for a symmetrical 3D null. They investigated the ef- 
fect of rotating the field lines around the spine and found that a 
twist wave (essentially an Alfven wave) propagates in towards 
the null. They found that while the helical Alfven wave spreads 
out, coupling due to the field geometry generates a fast wave 
that focuses on the null and wraps around it. This wrapping 
effect is analysed in more detail here for the simpler 2D null. 
The propagation of fast magnetoacoustic waves in an in- 
homogeneous coronal plasma has been investigated by 
[Nakariakov and Roberts (1995)[ who showed how the waves 
are refracted into regions of low Alfven speed. In the case of 
null points, it is the aim of this paper to see how this refraction 
proceeds when the Alfven speed actually drops to zero. 
The paper has the following outline. In Section |2] the basic 
equations are described. The results for an uncoupled fast mag- 
netoacoustic wave are presented in Section[3] Some simple an- 
alytical results are also discussed as a verification to and inter- 
pretation of the numerical simulation. Section H] discusses the 
propagation of Alfven waves and the conclusions are given in 
Section|5] 



where p is the mass density, v is the plasma velocity, B the 
magnetic induction (usually called the magnetic field), yu - 
An X lO^^Hm"' the magnetic permeability, 77 - liner is the 
magnetic diffusivity (m^s"'), and cr the electrical conductiv- 
ity. The gas pressure and the adiabatic energy equation are ne- 
glected in the low fi approximation. 

2.1. Basic equilibrium 

The basic magnetic field structure is taken as a simple 2D X- 
type neutral point. The aim of studying waves in a 2D config- 
uration is one of simplicity. The individual effects are much 
easier to identify when there is no coupling between the fast 
and Alfven modes. However, the extension to 3D is relatively 
straightforward. The modes will become coupled but their evo- 
lution is predictable from the 2D case. Therefore, the magnetic 
field is taken as 

Bo = Bo(-,0,--], (4) 

\a aJ 

where Bq is a characteristic field strength and a is the length 
scale for magnetic field variations. This magnetic field can be 
seen in Figure [T] Obviously this particular choice of magnetic 
field is only valid in the neighbourhood of the null point located 
at X = 0, z = 0. 
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Fig. 1. Our choice of equilibrium magnetic field. 



2. Basic Equations and Numerical IVIethod 



The usual MHD equations for a low j6 plasma appropriate to 
the solar corona are used. Hence, 



P\ 



dt 



+ (vV)v| = -(VxB)xB, 



dt 



5B , 
— = V X (v X B) -H T/V^B, 
dt 

(pv) = 0, 



(1) 
(2) 
(3) 



2.2. Linearised equations 

To study the nature of wave propagation near null points, the 
linearised MHD equations are used. Using subscripts of for 
equilibrium quantities and 1 for perturbed quantities, the lin- 
earised equation of motion becomes 



dt' 



VxB, 



xBo 



(5) 
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the hnearised induction equation 



3. Fast waves 



— i = V X (vi X Bo) + T/V^Bi , 
ot 

and the linearised equation of mass continuity 



dpi 
dt 



+ V-(poVi) = 



(6) 



(7) 



We will not discuss equation (Q further as it can be solved once 
we know vi . In fact, it has no influence on the momentum equa- 
tion (in the low /3 approximation) and so in effect the plasma is 
arbitrarily compressible (Craig & Watson, 1992). We assume 
the background gas density is uniform and label it as po. A spa- 
tial variation in po can cause phase mixing (Heyvaerts & Priest, 
1983). 

We now consider a change of scale to non-dimensionalise; let 
vi = vv;. Bo = BoB*), Bi = BoBp x = ax*, z = az\ V = iV* 
and f = If, where we let * denote a dimensionless quantity 
and V, Bq, a and t are constants with the dimensions of the 
variable they are scaling. We then set = v and v - alt 



(this sets V as a sort of constant background Alfven speed). This 
process non-dimensionalises equations Q and (|6]l, and under 
these scalings, f* = 1 (for example) refers to f = f = a/v; i.e. 
the (background) Alfven time taken to travel a distance a. For 
the rest of this paper, we drop the star indices; the fact that they 
are now non-dimensionalised is understood. 
The ideal linearised MHD equations naturally decouple into 
two equations for the fast MHD wave and the Alfven wave. 
The slow MHD wave is absent in this limit and there is no 
velocity component along the background magnetic field (as 
can be seen by taking the scalar product of equation (|5]l with Bo. 
The magnetic resistivity, tj, in equation (|6]l will be neglected in 
the numerical simulations but is included for discussion in the 
conclusions. 

The Hnearised equations for the fast magnetoacoustic wave are: 



dV 2, Jdb, db, 

— = V^ (X.Z) \ — — 



dt 
dbx 
dt 



vi(x,z) 

dV 
' dz 



\ dx 
db, 
'dt 



dz 
dV 
dx 



(8) 



where the Alfven speed, va (x,z), is equal to Vx^ + z^, Bi = 
(bx, 0, bz) and the variable V is related to the perpendicular ve- 
locity; V = [(vi X Bo) ■ Cyj. These equations can be combined 
to form a single wave equation: 



d^v 

5f2 



, Jd^V d^V 



(9) 



The linearised equations for the Alfven wave, with vi = 
(O, v,„ O) and Bj = (o, b,, o) are: 

dvy dby dby dby dvy dvy 
dt dx dz ' dt dx dz 
which can be combined to form a single wave equation: 



3h 
df 



d\ I d d 
'^'d-x-'d-zj 



(11) 



We solve the linearised MHD equations for the fast wave, 
namely equations (O, numerically using a two-step Lax- 
Wendroff scheme. The numerical scheme is run in a box with 
-6 < X < 6 and -6 < z < 6 and we initially consider a sin- 
gle wave pulse coming in from the top boundary. For the single 
wave pulse, the boundary conditions were set such that: 



V(x,6) 



sin for < f < - 
otherwise 



-^\jc=-6 - (J , -3-1. 

OX dx 



= 



dz 



Tests show that the central behaviour is largely unaflfected by 
these choices of side and bottom boundary conditions. The 
other boundary conditions on the perturbed magnetic field fol- 
low from the remaining equations and the solenodial condition, 
V-Bi =0. 

We find that the linear, fast magnetoacoustic wave travels to- 
wards the neighbourhood of the X-point and bends around it. 
Since the Alfven speed, va (x, z), is spatially varying, different 
parts of the wave travel at different speeds, and it travels faster 
the further it is away from the origin (i.e. the further away 
a point is from the origin, the greater in magnitude va (x, z) 
is). So the wave demonstrates refraction and this can be seen 
in Figure |2] A similar refraction phenomenon was found by 
INakariakov and Roberts (1995)| It is this refraction effect that 
wraps the wave around the null point and it is this that is the 
key feature of fast wave propagation. 

Since the Alfven speed drops to zero at the null point, the wave 
never reaches there, but the length scales (this can be thought 
of as the distance between the leading and trailing edges of the 
wave pulse) rapidly decrease, indicating that the current (and 
all other gradients) will increase. As a simple illustration, con- 
sider the wavefront as it propagates down the z axis at jc = 0. 
Here the vertical velocity is = ^ = -z- Thus, the start of 
the wave is located at a position - 6e^', when the wave is 
initally at z = 6. If the end of the wave leaves z = 6 at f = fo 
then the location of the end of the wave is Ze = 6e'°^' . Thus, 
the distance between the leading and trailing edge of the wave 
is 6z - 6{e'" - l)e ' and this decreases with time, suggesting 
that all gradients will increase exponentially. 



3.1. Current 

Since we have a changing perturbed magnetic field whose gra- 
dients are increasing in time, we have a build up of current 
density. Simulations show that there is a large current accumu- 
lation at the neutral point (Figure O and that this build up is 
exponential in time (Figure IH, in keeping with our discussion 
on the thickness of the wave pulse. Figure |4] shows the maxi- 
mum of the absolute value of the current at each time (left) and 
a plot of the logarithm of it against time (right). 
This build-up of current is extremely important since it implies 
that resistive dissipation will eventually become important, re- 
gardless of the size of the resistivity, and will convert the wave 
energy into heat. In fact, the exponential growth of the current 
means that the time for magnetic diffusion to become important 
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Fig. 2. Contours of V for a fast wave sent in from the upper boundary and its resultant propagation at times (a) t-l.O, (b) f = 1 .8, 
(c) t -2.6 and (d) f=3.4, labelling from top left to bottom right. Here, oj - In and time is measured in units of (jupo)'^^ /Bo- The 
axes have been rescaled in these lower two subfigures to draw attention to the central behaviour The number of points in each 
direction is 1600 x 1600. 



will depend only on log rj. Thus, refraction of the wave focuses 
the majority of the wave energy at the null point. This key re- 
sult will be discussed further in the conclusions. Note that the 
topology of the current accumulation seems to be approaching 
that of a current line. The current line comes from the collapse 
of the width of the fast wave as it approaches the null. Note 
that while the current grows exponentially in time, the velocity 
remains finite in magnitude. 



As an aside, note that the line in Figure|4](right) starts to deviate 
from a straight line after (approximately) t - 3.5. Recall the 
discussion above concerning the decrease in width of the pulse, 
where the distance between the leading and trailing edge of the 
wave was shown to be Sz = 6(e'° - l)e"'. This exponential 
decrease in length scales means that we will always run out 
of numerical resolution for the simulation in Figure |2l for any 
topology we want to resolve, for which 6z will be of the form 
^ (where is the number of grid points and C is a constant 
of inverse proportionality), we see that a solution is improperly 
resolved when e' > (e'» - 1)A/C,i.e. t > log A + log (e'° - 1)- 
log C. Thus, because of this logarithmic dependence on A, we 
will only be able to resolve the solution up until a certain time 
and increasing the number of grid points will not substantially 
improve this. Thus, the line in Figure|4]tails off due to a lack of 
numerical resolution. 



3.2. Analytical results 

We can approximately solve equation ^ for the fast wave to 
gain more insight into the current build-up observed in the nu- 
merical simulations. Substituting V - e'^*-*-^' ■ e^'"' into (|9]l and 
assuming that w » 1 (WKB approximation), leads to a first 
order PDF of the form 'F(x,z,(I>, - 0. Applying the 

method of characteristics, we generate the solution: 



z = 



X() cos 



zo cos 





. Axo 




+ zo sm 




\ Zo 




. (Axo 




- xo sm 


Zo I 


\ Zo 



-As 



where s is some parameter along the characteristic, xo is a start- 
ing point distinguishing between diff'erent characteristic curves, 
Zq is a second starting point (zo = 6 in our simulations), a> 
is the frequency of our wave and A is a constant such that 

A - Zo -^Xq + Zq. Figure |5] shows constant (p at four dif- 
ferent values of the parameter s. Constant can be thought of 
as defining the position of the leading edge of the wave pulse, 
i.e. with this choice of s, the WKB solution represents the front 
of the wave, s is comparable to t and so the subfigures can be 
directly compared to Figure|2l The agreement between the ana- 
lytic model and the leading edge of the wavefront is very good, 
as seen in an overplot of a contour and our WKB solution in 
Figure |6] Note that there is a diff'erence between the side and 
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Fig. 3. Shaded surfaces showing the build up of current at times (a) f=1.0, (fe) /=1.8, (c) f=2.6, id) f=3.4 and (e) t=A2, labelhng 
from top left to bottom right. Note the change in vertical scale. 




2 3 
tifne 



Fig. 4. ( - minimum current) against time elapsed (left), log ( - minimum current) against time elapsed (right). The slope of the 
Une between t= \ and ? = 3 is 1.65. 



bottom boundary conditions of the simulations and the analyt- from above can be approximated by: 
ical model. 



We can also use our WKB approximation to predict the current 
density build up. The current density is given by 7 = V which 
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Fig. 5. Plots of WKB solution for a wave sent in from the upper boundary and its resultant positions at times (a) t- 1 .0, (b)t- 1.8, 
(c) t = 2.6 and (d) t - 3.4, labelling from top left to bottom right. The lines represent the front, middle and back edges of the 
wave, where the pulse enters from the top of the box. Note the axes change in (c) and (d). 



By considering the modulus of this result, we can see that |j| 
will grow exponentially with s, i.e. as e^^\ and that the current 
will be negative (in agreement with figurelUi. To compare with 
our simulations, we first note that the time in our numerical 
scheme, t, is related to our parameter .? by f - cjs. Now, the 

exponent 2As in the exponential is equal to llusj + and 
so along jcq = and at zo = 6 (xq = is where the maximum 
current of our numerical simulations occured) this is simply 2f . 
The slope of the numerical experiment in figure|4]between t - I 
and f = 3 is 1.65. This agreement between the analytical and 
numerical current density build up is quite good, considering 
that the WKB solution is only valid for a harmonic wavetrain 
with w » 1 . 

3.3. Asymmetric Pulse 

As mentioned in the introduction, MHD waves in the neigh- 
bourhood of null points have been investigated before in cylin- 
drical coordinates. In this paper however, we can readily con- 
sider asymmetric pulses and how they propagate in the neigh- 
bourhood of coronal null points. Hence, to demonstrate the 
wrapping around the null phenomena further, we consider the 
same numerical experiment above but with a pulse being fed 
in only on one part of the upper boundary. This pulse took the 
form; 

y = sin(o.0[l+cos7r(x-3)] |o^^|? '^"^ (12) 



The propagation can be seen in Figure|2l Again, it is important 
to note the majority of the wave energy tends to accumulate at 
the coronal null point; this is where wave heating would occur. 
The growth of the current density occurs in a similar manner to 
that demonstrated in Figure[3] However, it is important to note 
that the more the pulse is displaced to one side, the more the 
boundary conditions will influence the subsequent evolution. 

4. Alfven waves 

The equations describing the behaviour of the Alfven wave, 
equations ( fTOl i. were solved numerically using the same two- 
step Lax-Wendroff scheme. We initially consider a box (0 < 
X < 6, < z < 6) with a single wave pulse coming in across 
half of the top boundary (0 < x < 3). We chose such a pulse be- 
cause, as shown in Figure[8l the Alfven wave spreads out along 
the field lines as it propagates and we found that this choice of 
boundary condition illustrated this effect much clearer. The full 
boundary conditions were; 

/ • irx\ I, f < X < 3 

Vy(x,o) = sm wr(^l -I- cos -J- J for i q < ^ < ;r 
-^\z=6 = otherwise 

dVy dVy dVy 

^l.v=6 = , -^|,=o = , -^|,=o = . 

ox ox oz 

Tests show that the central behaviour is unaffected by these 
choices. The other boundary conditions follow from the re- 
maining equations and the solenodial condition, V ■ Bi =0. 
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Fig. 6. Comparison of numerical simulation and analytical solution at times (a) f=1.0, {b) t - 1.8, (c) t - 2.6 and {d) t - 3.4, 
labelling from top left to bottom right. As before, the lines are the front, middle and back edges of the analytical solution. 



Note that we have used a slightly different inital pulse to those 
in the fast wave investigation. This is because the field lines 
(see Figure[Tl) leave the box and we know that the Alfven wave 
follows the field lines. Hence, we made this choice of inital 
pulse in case reflections from the side boundaries influenced 
the subsequent evolution. 

We found that the linear Alfven wave travels down from the top 
boundary and begins to spread out, following the field lines. 
As the wave approaches the lower boundary (the separatrix), it 
thins but keeps its original ampitude. The wave eventually ac- 
cumulates very near the separatrix; defined by the x axis. This 
can be seen in Figure[8] 

4.1. Analytical results 

The Alfven equation we have to solve takes the form: 

d\ 



dx ^ dzj 



and this can be solved using the method of characteristics. Let 
Ts ~ {■"■Jk ~ ^i^) '^^'^ comparing the original equation with 



^ = ^ ^ + 1^ % leads to: 

OS OS ax OS oz 



Xi)e 



Z{)e 



(13) 



where xq and zq are the starting positions of our characteristics. 
In our simulation, zq - 6. Thus, our characteristic equation, 

-Tjr- = can be solved with a D' Alembert solution such that: 

oT-^ as- 



In order to compare these analytical results with the numerical 
results above, we substitute the same initial conditions into the 
D' Alembert solution, i.e. T (xq) - 1 + cos(^) and g(t) - 
sin (ojt) to get the analytical solution for Vy, namely: 



Vy(X, Z, t) = 



1 + COS 



for 



(TTXZW . ( , Z\ 

sm w f + log - 

I 18 /J \ ^6/ 



< f + log § < i 
0< 2^ < 3 



(15) 



The agreement between the analytical and numerical results 
is excellent (the contours essentially lie on top of each other), 
even though the analytical solution does not satisfy the numer- 
ical boundary conditions. 

Furthermore, we use our analytical solution to calculate by, 
and ;,: 



by = 

jx = 
jz = 



2k 

z 



1 +cos(^j sinw^f + log|j , (16) 
(^|)]cos^(r + log|) 

n(^)sinw(f -i-log|j , (17) 

jsin w^f -I- log |j (18) 



1 + cos 



nx . 

sm 

18 



7TZ . jnxz 

— sm 

18 I 18 



Vy^T(xQ)g{t-s) . 



(14) 



all for 



0<f + log|<S 
< ^ < 3 
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Fig. 7. Contours of V for a fast wave sent in from upper boundary in accordance with equation ( fT2b with o) -2n , and its resultant 
propagation at times (a) t= 0.5, {b) t= 1.17, (c) t = 1.83, {d) t- 2.33, (e) f= 3.17 and (/) t- 3.67, labelling from top left to bottom 
right. The axes have been scaled down in the last two subfigures to draw attention to the central behaviour. The number of grid 
points in each direction is 1200 x 1200. 



4.2. Current 

As in the fast wave case, we have a spatially varying perturbed 
magentic field and so current is forming, given by ^ (V x Bi) = 

{jx, 0, jy). In the Alfven case only - ~j;''j7 ™d ~ jlTf 
present. The evolution of the current can be seen in Figures |9] 
andim 

From Figure |9l we see that spreads out along the field lines, 
accumulating along the separatrix; z = 0. /i takes a discontin- 
uous form; this is due to our choice of initial conditions and 
this is confirmed by the analytical solution, jx also increases in 
time. From equation[T7] we see that jx grows like i and thus, in 
accordance with equation [13] this means jx grows like e', i.e. 
grows as e' (since s - t + constant in the D' Alembert solution). 
The behaviour of the maximum of jx with time can be seen in 
Figure [TOl the slope of the line is + 1 .0 in agreement with our 
analytical solution. 

The behaviour of j, can be seen in Figure [TT] takes the spa- 
tial form of the inital pulse and, in the same way as y,., spreads 
out along the fieldlines. It decays in amplitude as it approaches 
the separatrix. The analytical solution, equation[T8] shows that 
it behaves like z (and z is going to zero). Thus, according to 
equation [13] j, decays as e^'. The behaviour of the maximum 
of with time can be seen in Figure[T2] In this case, the slope 
of the line is -1.0 in agreement (again) with our analytical so- 
lution. 



Hence, the Alfven wave causes current density to build up 
along the separatrix. 

5. Conclusions 

This paper describes the start of an investigation into the nature 
of MHD waves in the neighbourhood of null points. From the 
work explained above, it has been seen that when a fast mag- 
netoacoustic wave propagates near a magnetic X-type neutral 
point, the wave wraps itself around the null point due to re- 
fraction (at least in two dimensions). It has also been seen that 
this behaviour causes a large current density to accumulate at 
the null and simulations have shown that this build up is expo- 
nential in time, although the exponential growth in this linear 
simulation will be modified by non-linearities. We also note 
that for the set of disturbances investigated here, there is no 
evidence of the X-point collapsing; rather, the current density 
seems to form a spike. However, it is clear that the refraction 
of the wave focusses the energy of the incident wave towards 
the null point. As seen from both the numerical work and ana- 
lytical approximations, the wave continues to wrap around the 
null point, again and again. The physical significance of this 
is that any fast magnetoacoustic disturbance in the neighbour- 
hood of a neutral point will be drawn towards the region of zero 
magnetic field strength and focus all of its energy at this point. 
Hence, this is where the build up of current will occur and en- 
ergy will be dissipated. Experiments are being carried out to ex- 
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Fig. 8. Contours of Vy for an Alfven wave sent in from upper boundary for < x < 3.0 and its resultant propagation at times (a) 
f=0.25, {b) f=0.5, (c) t=Q.15, (d) f =1.0, (e) t=1.25 and (/) t=1.5, labelling from top left to bottom right. 



tend the analysis to multiple null points and three dimensions 
but, if the results transfer, then null points should effectively 
trap and dissipate the energy contained in fast magnetoacoustic 
waves. Therefore, wave heating will naturally occur at coronal 
null points. 

The numerical experiments and analytical work described 
above were all conducted using the ideal MHD equations. 
However, we can make some comments about the addition 
of resistivity into the model. For the fast magnetoacoustic 
wave, all the current density accumulates at the null point 
and appears to form a null line. Hence, no matter how small 
the value of the resistivity is, if we include the dissipative 
term, then eventually the t^V^Bi term in equation (|6]l will 
become non-negligible and dissipation will become impor- 
tant. In addition, since V-Bi grows exponentially in time, 
the diffusion terms become important in a time that de- 



pends on log 77; as found by Craig and Watson (1992) and 



Craig and McClymont (1993) This means that linear wave dis- 
sipation will be very efficient. Thus, we deduce that null points 
will be the locations of wave energy deposition and preferential 
heating. 

In the case of the Alfven wave, the results show that the wave 
propagates along the field lines, accumulating on the separatrix 
and hence, due to symmetry, along the separatrices. The wave 
also thins and stretches along the separatrices. The current 
increases and accumulates along the separatrix, whilst j- de- 
cays away. This is seen in both the analytical and numerical 
work. 



Now consider the effect of including resistivity. Considering 
V^B 1 , this can be split into # ^ -H # ^ = |^ -H ^ . As shown in 

^ ' ^ ox ox oz oz ox oz 

the appendix, ^ decays away (exponentially) but ^ increases 
(exponentially). We have also seen that the current accumulates 
along the separatrices. Hence, eventually (due to its exponential 
increase), the resistive term, r]'-^ - rj'^, will become impor- 
tant, no matter how small the value of 77. Hence, all the Alfven 
wave energy will be dissipated along the separatrices. This is a 
different behaviour to that of the fast wave in the sense that the 
two wave types deposit all their wave energy at different areas 
(along the separatrices as opposed to the null point), although 
the phenomenon of depositing wave energy in a specific area is 
common to both. 

One of the next steps to be taken is to extend the model to 3D. 
Galsgaard et al. (2003) performed a similar analysis of MHD 
wave propagation for a particular disturbance of a symmetric 
3D null. Their findings show a close link to our own and sup- 
port the possibility of transferring our results to 3D. 

Another of the steps to be taken is to investigate the ef- 
fect of pressure on the system. This has been investigated by 
Craig and Watson (1992)[ with cylindrical symmetry and they 
find that the rapid current growth is halted. The most obvious 
effect of including a finite fi is the introduction of slow mag- 
netoacoustic waves. Fast waves can now pass through the null 
point (as we would now have a non-zero sound speed) and thus 
perhaps take wave energy away from that area. The exact na- 
ture depends on the choice of boundary conditions. In a simple 
manner, if the yS has a value of ySo ^ 1 at, say, + - 1, the 
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finite pressure effects will become important once - 2/3. 

From our WKB solution, with zo - 6 and xq - 0, the gas 
pressure will become non-negligible when t - - \ logjSo. This 
is not true for Alfven waves since they are unaffected by the 
finite y6. 

Finally, the validity of the linearisation is questionable once 
the perturbed velocity becomes comparable to the magnitude 
of the local Alfven speed. In a similar argument to the finite fi 
case, the linearisation is valid until t ^ -logM^i, where Ma is 
the initial Alfven Mach number. Once the Alfven Mach num- 
ber exceeds unity, the fast wave is likely to shock and the wave 
energy will still be dissipated but this time in the shock. A more 
rigorous analysis requires detailed numerical simulations of the 
non-linear MHD equations and will be investigated in the fu- 
ture. 
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Appendix 

Consider the induction equation for the Alfven wave: 
dby Id d\ d d 

In order for the resistive terms to become important, they have 
to become comparable to the convective term; [x-^ - Vy. 
So we have to calculate the behaviour of j^j^ and the behaviour 
of yJx over time. From Figure[T3l we see that the maximum of 



j, decays exponentially over time. From Figure [141 we see 
that the maximum of -4;]^ increase exponentially over time. 



d: 

dx 

Hence, the resistive terms invloving -^j^ will eventually be- 
come important, no matter how small the value of the resis- 
tivity, whereas the resistive terms invloving -^j^ will become 
negligible. 
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Fig. 10. Maximum jx against time elapsed (left), log ( maximum jx) against time elapsed (right). The slope of the line is +1.0. 
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Fig. 12. Maximum against time elapsed (left), log ( maximum j^,) against time elapsed (right). The slope of the line is -1.0. 
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Fig. 13. Plot of maximum 4-]^ against time (left), plot of logmaximum^jj. against time (right). The slope of the hne is -2 




